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Abstract 

By means of the Furnstahl, Serot and Tang's model, the effects of surface 
tension and Coulomb interaction on the liquid-gas phase transition for finite 
O ! nuclei are investigated. A limit pressure p\\ m above which the liquid-gas phase 

transition cannot take place has been found. It is found that comparing to 
the Coulomb interaction, the contribution of surface tension is dominate in 
low temperature regions. The binodal surface is also addressed. 
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Since the arguments given by Miiller and Serot [1] that a second order liquid-gas (L-G) 
phase transition will take place in a multi-components and multi-conserved charged system, 
much theoretical efford has been devoted for studying this problem by using different models 
and different treatments [2-4]. But all investigations are limited to infinite nuclear matter. 
It is of interest to extend this study to finite nuclei. This is the objective of this paper. 

If we consider the finite nuclei as a liquid droplet and discuss its L-G phase transition, 
two major effects, namely, the surface energy of the droplet and the Coulomb interaction 
of proton-proton must be considered. The reasons are as follows: it has been shown that 
the difference of the chemical potentials between proton and neutron plays an essential role 
to determine the order of L-G phase transition [I, 2]. For the infinite nuclear matter the 
chemical potentials of proton and neutron depend on the third component /3 of isospin when 
nucleon-nucleon-p-meson (NNp) interaction exists. In symmetric nuclear matter, the L-G 
phase transition is of first order because I3 = 0. In asymmetric nuclear matter, 1 3 ^ and 
then the chemical potential of neutron p n does not equal to the chemical potential of proton 
fi p , a second order phase transition may take place. The Coulomb interaction cannot be 
taken into account because it becomes divergent in infinite nuclear matter. But in finite 
nuclei, the contribution of Coulomb interaction can be considered. Obviously, the chemical 
potential of proton p p not only depends on J 3 , but also on Coulomb interaction. But the 
later has no effect on fi n . The contribution of Coulomb interaction will make that the values 
of fip and fi n becomes more different. 

Besides Coulomb interaction, on the other hand, the surface tension of the droplet will 
affect on the pressure of the liquid phase and then on the coexistence equations because the 
pressures of two phases must equal at the phase transition point. 

To exhibit the effects of surface energy and Coulomb interaction on the L-G phase tran- 
sition for finite nuclei, we employ the Furnstahl-Serot-Tang (FST) model [5-8], which has 
been shown to be successful to explain the properties of both infinite nuclear matter and 
finite nuclei. The Lagrangian density of FST model under mean field approximations reads 
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where g s , g v g p are, respectively, the couplings of light scalar meson a, vector meson u and 
isovector meson p fields to the nucleon, <f> , V , bo are the expectation values 0o =< >, 
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The scalar fluctuation field <p is related to S by S (x) 



So — 4>{x) and H q is given by m s = AH q / (d Sq), d the scalar dimension. By using the 
standard technique of statistical mechanics, we get the thermodynamic potential Q as [2, 9] 
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where f3 = 1/ksT and the quantity z/j (i = n,p) is related to the usual chemical potential pi 
by the equations 
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where p 3 = p p —p n , the third component of isospin J 3 = (N p — JV n ) /2 = Vp 3 /2, and _E* (fc) = 
\JM* 2 + /c 2 with M* = M — g s 0o- Usually, instead of p 3 , we introduce ct = (p n — p p ) / p, the 
asymmetric parameter to calculate, where p = p n + p p . Having obtained the thermodynamic 
potential, all other thermodynamic quantities, for example, pressure p, can be calculated. 
We get 
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are the nucleon and anti-nucleon distributions respectively. The neutron density p n and the 
proton density p p are given by 
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Now we are in a position to study the L-G phase transition of finite nuclei. The two-phase 
coexistence equations are 
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where the prime and the double prime refer to the liquid phase and gas phase, respectively. 
Considering the droplet as an uniformly charged sphere, the contribution of Coulomb inter- 
action to the chemical potential of proton and the pressure are [9, 10, 11]. 
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respectively. The additional pressure provided by the surface tension of the liquid droplet 

is [12, 13] 



Psur f (T,p) = -2 1 (T)/R (14) 
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with T c being the critical temperature of L-G phase transition in symmetric nuclear matter. 
We take the liquid droplet along the (3 -stability line, it satisfies 

Z = 0.5A - 0.3 x 10~ 2 A 5/3 (16) 

The parameters be chosen for our numerical calculations are the set T\ of FST model 

^ = 99.3, tf = 154.5, ^ = 70.2 (17) 

m s = 509MeV, S = 90.6MeV 
C = 0.0402, r] = -0.496, d = 2.70 

Our results are summerized in Fig.l-Fig.6. To make our results more transparent, we 
neglect the Coulomb interaction and the surface effect by taking R — > oo in Fig.l and the 
p n , Up isobar vs a reduce to that of infinite asymmetric matter. In this Figure, we fix the 
temperature T = 10 MeV and the curves a, b, c, d, e correspond to the pressure 0.06, 0.085, 
0.100, 0.164 and 0.200 MeV/m" 3 respectively. We see that the curves for lower pressures 
are more complicate than those of the large pressures. When p = 0.200 MeV/m~ 3 , curve 
e has one branch only, but when p = 0.06 MeV/m" 3 , curve a has three branchs. The 
chemical potentials isobar /i", \j!' vs a" given by the right hand side of eq.(9-ll) for the 
gas phase is shown in Fig. 2 where T = 5 MeV and p = 0.016 MeV/m -3 . In fact, these 
curves are the same as that of infinite nuclear matter because the chemical potential and 
the pressure for the gas phase do not depend on the Coulomb interaction and the surface 
term. We see in this case both pi' (a") and p!' n (a") curves have three branchs. The chemical 
isobar as a function of a for the liquid phase and the gas phase are shown in Fig. 3 by solid 
line and dashed line respectively where we fixed T = 5 MeV and p = 0.016 MeV/m" 3 . 
The dashed lines in Fig. 3 for gas phase are in fact the same curves as those of Fig. 2 except 
that the range of the ct-axis is (0.0, 0.5) instead of (0.0, 1.0). The rectangle construction [1, 
2] which represents for the Gibbs' conditions of eq.(9 - 11) for the two-phase equilibrium is 



also plotted in Fig.3 by the edges of a rectangle. Due to the effect of Coulomb interaction 
and the surface energy, the chemical potential isobars for the gas phase and for the liquid 
phase are very different. We see from Fig.3 that the p," (a") and p!' n (a") curves for the gas 
phase have two branchs shown by dashed lines in the regions < a < 0.5 respectively, but 
p' (a') p' n (a') for the liquid phase has one branch in this region only. This behavior is quite 
different from that of infinite nuclear matter in which the liquid phase and the gas phase 
chemical isobars p n and p p are shown by the same curves (for example , see curve a of 
Fig.l). This new feature leads to multi-solutions for Gibbs' condition because one can find 
four rectangles between two branchs of p!' n ip'p) and one branch of p' n ( p' p ) . But according 
to the equilibrium condition: the chemical potential of the system in equilibrium state must 
take the minimal value at fixed temperature and pressure [14]. Therefore, only one rectangle 
which corresponds to the minimum chemical potential isobar shown by solid lines in Fig.3 
refers to a stable equilibrium phase transition, and the others are all metastable states. The 
other three rectangles are not shown in Fig.3. 

The section of binodal surface [1,2] at finite temperature T = 5 MeV is shown in Fig. 4. 
A limit pressure pn m = 0.018MeV/m~ 3 above which the rectangle cannot be found and the 
coexistenced equations eq.(9-ll) have no solution has been obtianed. The binodal surface 
will cut off at limit pressure pn m . This situation is very similar to that of our previous paper 
[2] in which we considered the density dependence of the NNp coupling g p (p). The reason 
is that no matter g p (p) or the Coulomb interaction or the surface energy, even though 
they change the chemical potential \x n and \x v in different fashions, they will make that 
the rectangle construction turns out to be disappear. In fact, this result is a reflaction of 
the so called Coulomb instability in finite nuclei [9-11, 13, 15]. The Coulomb instability 
of FST model has been discussed in detail by our previous paper [9] in which we found a 
limit temperature T\ ira above which the coexistenced equations have no solution and the 
L-G phase transition can not take place. The difference is that, instead of T\ ira , we now fix 
temperature T = 5 MeV to find the limit pressure. 

Finally, we hope to compare the effects of Coulomb interaction and the surface tension on 
the L-G phase transition, separably. As shown in eq.(13) and eq.(14), we see that, firstly, the 
effects of Coulomb interaction and the surface tension are opposite because of pcoui and p sur f 
with opppsite signs, and secondly, p sur f depends on temperature but pcoui is independent. It 
means that the Coulomb interaction and the surface tension play different roles in different 
temperature regions. In low temperature regions T<5 MeV<^cT c , 7 (T) becomes larger, and 
we have |p sur /| > Pcouh the surface tension becomes dominant. To show this result clearly, 
we draw the chemical isobar vs. a curves with the surface effect only and with the Coulomb 
interaction only in Fig.5 and Fig. 6 respectively. Comparing Fig. 5, Fig. 6 and Fig.3 we see 
that the curves of Fig.5 are very similar to those of Fig.3, but the curves of Fig.6 are very 
different. It confirms the surface effect dominates at T = 5MeV. The curves in Fig.6 are 
very complicate. The reason is that we see from eq.(ll), in the phase transition process, 
Psurf "will increase p" since it has negative sign, but p' Cou i will decrease p" . As indicated by 
Fig.l, the pressure p" becomes lower, the curve becomes more complicate. 

In summary, it is shown that the surface effect and the Coulomb interaction are important 
for the L-G phase transition of finite nuclei. In low temperature T <^T C regions, the surface 
effect is dominate. A limit pressure pu m above which the L-G phase transition cannot 
take place has been found. Since the critical temperature T c = 14.75 MeV, and the limit 



temperature Tn m is around 5.4 — 8.8 MeV for FST model [12], we come to a conclusion that 
the surface effect is dominate in the L-G phase transition of finite nuclei for FST model. 
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FIGURES 

FIG. 1. The chemical isobars for infinite nuclear matter, where T = 10 MeV, and a, b, c, d 
e refer to the pressure 0.06, 0.085, 0.100, 0.164 and 0.200MeV/m~ 3 respectively. 
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FIG. 2. The chemical isobars for the gas phase where T = 5 MeV and p = 0.016 MeV/m 3 . 

FIG. 3. The rectangle construction for two-phase equilibrium for T = 5 MeV and p = 0.016 

MeV/m" 3 . 

FIG. 4. The section of binodal surface for T = 5 MeV. 

FIG. 5. The chemical isobars vs. a curves with surface effect only. 

FIG. 6. The chemical isobars vs. a curves with Coulomb interaction only. 
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